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DETERMINATION OF ELASTIC STRESSES IN GAS-TURBINE DISKS 

By S. S. Manson 


SUMMARY 

A method is presented for the calculation of elastic stresses 
in symmetrical disks typical of those of a high-temperature gas 
turbine. The method is essentially a finite-difference solution 
of the equilibrium and compatibility equations for elastic 
stresses in a symmetrical disk . Account can be taken of point- 
to-point variations in disk thickness, in temperature, in elastic 
modulus, in coefficient of thermal expansion, in material density, 
and in Poisson’s ratio. No numerical integration or trial-and- 
error procedures are involved and the computations can be per- 
formed in rapid and routine fashion by nontechnical computers 
with little engineering supervision. Checks on problems for 
which exact mathematical solutions are known indicate that the 
method yields results of high accuracy. 

Illustrative examples are presented to show the manner of 
treating solid disks, disks with central holes, and disks con- 
structed either of a single material or of two or more welded 
materials. The effect of shrink fitting is taken into account by 
a very simple device. 

INTRODUCTION 

One of the problems in the design of gas turbines is the 
determination of the stresses in the turbine disk under 
operating conditions. Calculation of the elastic-stress dis- 
tribution is a first step in the determination of the true stress 
distribution. This stress distribution is based on the assump- 
tion of linearity of stress with strain and differs from the true 
stress distributions, which may contain stresses beyond the. 
proportional elastic limit of the material. 

The equations for the elastic-stress distribution in sym- 
metrical disks are well known. Their solution may, however, 
offer considerable difficulty. One difficulty encountered in 
calculation of the operating stresses in disks with liigh- 
temperature gradients is that the physical properties of the 
materials, such as elastic modulus, Poisson's ratio, and 
coefficient of thermal cipansion, vary with the temperature 
and therefore have a different value at each location in the 
disk. In addition, the thickness of the turbine disk varies 
from radius to radius. If the disk consists of portions of 
different materials welded to each other, the density may 
vary from one section to the other. Shrink fitting and weld- 
„ ing of the component parts at elevated temperatures also 
introduce special stress problems. Attempts to find complete 
analytical solutions for the stress problems that take into 


account shrink fitting and point-to-point variation of the 
physical properties and of the disk thickness result in math- 
ematical complexities; approximate solutions are therefore 
usually found by numerical methods. 

Thompson in reference 1 gives a numerical approach to the 
turbine-disk problem that takes into account point-to-point 
variation in disk thickness and in all physical properties 
except Poisson’s ratio. A method capable of easily account- 
ing for shrink fitting and for the variation in Poisson’s ratio 
as well as in the other properties was developed in 1945 at 
the NACA Cleveland laboratory and is presented in the 
analytical section herein. The method is essentially a finite- 
difference solution of the differential equations of stress in a 
rotating disk and incorporates several advantageous features 
uncommon to other forms of solution. For example, numer- 
ical integration and trial-and-error processes have been 
completely avoided, which makes it possible for nontechnical 
computers to carry through the entire solution rapidly and 
with little engineering supervision. 

In the second section of this report, illustrative examples 
are presented to show the manner of treating a solid disk and 
a disk with a central hole for application to a gas turbine. 
The examples are self-explanatory and may be used as a 
guide for the solution of turbine-disk problems without 
reference to the analytical section of the report in which the 
basis for the solution is derived. 

The method also has applicability to the study of stresses 
in rotating disks other than that in the gas turbine. It has 
been applied, for example, to disks of axial-flow compressors 
in which the only complicating factor is variable disk thick- 
ness. For such application the main advantages of the 
method are the routine nature of the. calculation, the rapidity 
with which the calculations can be made, and the accuracy 
of the final results. 

ANALYSIS 

SYMBOLS 

The following symbols are used: 

A point midway between nth and (n-l)st point stations 
E elastic modulus of disk material, (lb/sq in.) 
h axial thickness of disk, (in.) 

r radial distance, (in.) 

u radial displacement of any point on disk as disk 
passes from unstressed to stressed condition 
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a coefficient of thermal expansion between actual tem- 
perature and temperature at which there is zero 
thermal stress, (in./(in.)(°F)) 

A T temperature increment above that at which there is 
zero thermal stress 
6 r radial strain, (in. /in.) 

t t tangential strain, (in./in.) 

At Poisson’s ratio 

p mass density of disk material, ((lb)(sec 2 )/in. 4 ) 

(t t radial stress, (lb/sq in.) 

( r t tangential stress, (lb/sq in.) 

co angular velocity of disk, (radians/sec) 

The following supplementary subscripts are used for 
denoting values of the preceding symbols in connection with 
the finite-difference solution: 
n nth point station 

7i-l (n-l)st point station 

a station at the smallest radius of the disk considered 

(For a disk with a central hole, this station is taken 

at the radius of the central hole ; for a solid disk, 
this station is taken at a radius of about 5 percent of 
rim radius.) 

b rim of disk or base of blades 

Example of the use of double subscript: 
ar.n-i radial stress a T at station n—\ 

The following supplementary symbols denote combina- 
tions of the foregoing symbols arising in the analysis: 

-A T ,n 


A. t,n 
B T ,n 

Bt.n 
C V nb/n 


coefficients defined by equations 
(Tr.n — ^lr.n (Te,a~i~ B T<n 

(T B t,n 


ru 

C n ~E n 


( 1 ~f~ Pti) iTn ^n —l) 
2Er,r n 


E n cy (r n r n _j)A n 


D' n 


1 i (1 T Mn) (^n i) 


E n 1 

F « — Tn-xh-n—i 


2E n r n 


jp/ __ Pn—1 (l~f" Mn — l) (r n T n - 1 ) 

n ~E»- 1 " ' 2E n _ir n . i 


G n 2 (? n ^n—\)btn—l 

l (l + M,-l)(r B 


G\= 


■r n - 1 ) 


E n — i 2E n -\T n _\ 

1 


Fin 2 03 G n ^n—\) (^PuF-nTn Pn— \^n—Fn—\^ 
H n ' = (Xn^T n OC n ~i&T n _i 

F' T) —FB' 

tt- * n 1 -' n x n 1 -' n 

— pi t\ _ ri jy 

' y n^n n 


K' n 


L n =~ 


C n F' n — C' n F n 

C'Jh-CnD'n 

G' n D n +G n D' n 
C' n D n —C n D' n 


L'n- 


M n 


C' n G n +C n G' 

C'„D n -C n D' 

H' n D n +H n D' n 

C' n D n -C n D' n 


n 

n 


i// C' n H n -\-C n H' n 
n ~C\D n -C n D' n 


ASSUMPTIONS 

The assumptions are made that stress is proportional to 
strain and that the disk material is completely elastic at the 
stress distribution induced by the centrifugal and thermal 
effects. All variables of material properties and operating 
conditions are assumed to be symmetrical about the axis of 
rotation. Axial stresses are neglected and at any radius the 
radial and tangential stresses are assumed to be uniform 
across the thickness of the disk. Temperatures are taken in 
the central plane perpendicular to the axis of the disk. 

OUTLINE OF METHOD 

In a thin rotating disk of variable thickness, the state of 
stress at any radius can be completely defined by the two 
principal stresses, the radial and tangential stresses cr r and a t , 
respectively. Two equations are therefore necessary to deter- 
mine the two unknown stresses. The first of these equations 
can be obtained from the conditions of equilibrium of an 
element of the disk; the second, from the compatibility con- 
ditions, which are mathematical statements of the interrela- 
tion between the radial and tangential strains in a sym- 
metrical disk. 

The equilibrium and compatibility equations result in 
differential form defining relations between the stresses at 
radius r and those at a radius infinitesimally removed from r. 
Except for some special cases, the solutions of these equations 
are difficult to obtain. In order to facilitate solution, the 
differential equations are rewritten in finite-difference form 
relating the stresses at radius r with those at a radius finitely 
removed from r. By means of the finite-difference equations, 
the stresses at an arbitrary finite number of stations along 
the disk radius are expressed in terms of the stresses at a 
single reference station near the center of the disk. For a 
disk with a central hole the reference station is chosen at the 
inside radius, where the radial stress is zero; hence, the 
stresses at all stations in the disk are expressed in terms of the 
single unknown, the tangential stress at this station. For a 
solid disk, the reference station is chosen at a point near the 
center of the disk (at a radius of about 5 percent of the disk 
radius). In this region the radial and tangential stresses can 
be assumed to be approximately equal; again, therefore, the 
stresses at all stations are expressed in terms of a single 
unknown. The unknown can then be determined by the 
boundary conditions at the rim of the disk where the radial 
stress is equal to the centrifugal blade loading. When the 
radial stress at the rim, expressed in terms of the tangential 
stress at the reference station, is equated to the blade loading, 
the tangential stress at the reference station is evaluated. 
After the tangential stress at the reference station has been 
determined, all the other stresses, expressed in terms of this 
stress, can be evaluated. 
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DIFFERENTIAL EQUATIONS 


The equilibrium equation, as given in reference 2 (p. 374), 
using the notation of this paper, is 

£.(rh<,,)-h<r,+l>Mr*=0 (1) 


The compatibility relation is obtained by elimination of u 
from the stress-strain displacement equations 


e r = 


du 

dr 



-\-a.AT 


( 2 ) 


t t =~-~-£° r + aAT (3) 


Equation (3) is subtracted from equation (2) to eliminate 



du u (1 T m) (<r r— (T t) 
dr r E 

(4) 

or 

rdu~udr (l-}-j u)(o>— <r £ ) 
rdr E 

(5) 

But 

rdu—udr d /w\ 
r' 2 dr dr \r) 

(6) 


Therefore, by equations (3), (5), and (6), 


d_ ( — 

dr\E) dr \E ) 


i(^T) 


(1+ fl)(<T r ~ <rQ _ n 

Er " U 


(7) 


Equations (1) and (7), together with a knowledge of the 
boundary conditions, are sufficient to solve for the two 
unknowns a r and a,. Because p, E, p, a, A T, and h are, in 
general, functions of the radius r, the equations cannot 
readily be solved in their differential form; a finite-difference 
solution was therefore derived. 


FINITE- DIFFERENCE EQUATIONS 

The translation of differential equations into finite- 
difference form to facilitate solution is common in engineering 
practice. The method has, in fact, been applied in limited 
fashion to the solution of the steam-turbine disk problem 
(reference 2, pp. 398-400). This application neglects, how- 
ever, the point-to-point variation in physical properties, and 
therefore no application to the gas-turbine disk, in which 
there is appreciable variation in properties from hub to rim, 
is made. In addition, the solution of the equations involves 
an inte”- dation procedure, the elimination of which could 
considerably reduce the amount of calculation necessary for 
a solution and increase the accuracy of the final results. 

A number of discrete point stations are chosen along the 
disk radius as shown in figure 1 (a). If it is assumed that the 
stress distribution in the disk has already been determined, 
all quantities appearing in equations (1) and (7) are known 
at each of the point stations and the values of corresponding 
quantities at the point A midway between the nth and (n-l)st 
point stations can be approximately determined. For 



(a) Location of point stations. 

(!>) Value of typical function midway between point stations n and n- 1. 

Figure 1. — Sketches used to derive finite-difference equations for stresses in symmetrical 

rotating disk. 


example, in the plot of rha T against r (fig. 1 (b)), the radius at 
point A is expressed as 

r A = \ (r„-i + r n ) 

the value of rha T is 

(^rfl<rd) A ~ (^n — l^n— l® - r,n— 1 “f” r n^n^r.n) 


and the slope of the curve at point A, which is approximately 
equal to the slope of the chord joining points n and n- 1, is 
defined as 


d / i \ r n }l n (TT,n ) n—\h'n—\ ( r r,n — l 

dd rU; >*“ 


In a similar way, the values of each of the other variables 
entering into the equations can be evaluated for point A. 
If the evaluations are correct, the quantities at A must satisfy 
equations (1) and (7). These equations therefore become, in 
finite-difference form 

r nh'n.Gr.n lOV.n— 1 hn& t, n ~f~ ^ n — 1 & t,n — 1 _|_ 

r n —r n - 1 2 

2 (.Pnh'nt'n~ T Pn— ill n — U n — l”) 0 (8) 
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and 

G t, n G l,n — 1 n 1 G r ,n— 1 

/?n-i _ E n E n _x , a n AT„-a B _ 1 A !T B . 1 _ 

f’n—rn-i r n —r n -i r n —r n _ x 

1 [~ ( 1 4~M«) (°V,ra G t,n) , (1 ~f" Mrc— l) (o"r,«— 1 & t,n — 1)~ | A 

2[_ E n r n + W 

which reduce to 

C n<S r,n B n <7 t,n~ F nGr,n—l~\~ G n (T t,n — 1 Hn (10) 

and 

C n<f T . n — D' n a t, n =F r n a Ttn _ x — G' n cr t , n -iFH' n (11) 


of elasticity (1) and (7), or their ultimate finite-difference 
forms in equations (13) and (14). If an equation in the form 
cx + d — 0 is to be true for all values of x, the coefficients 
c and d both must be zero. Because equations (13) and (14) 
are to be true, independent of the value of cr,, a , the coeffi- 
cients of cr, „ must be zero, and the two equations reduce to 

GjiA r , n F) n At, n F nAf, n— 1 G n At n — 1 — 0 \ 

C' n A T , n — D' n A t , n — F' n A r , re _i + G' n A t , n _i — 0 1 

/ (15) 

C n B r , n — D n B t , n — F n B r „_i — G n B t . „..i-i-^ re =ol 
C' n B r , n ~D' n B t .n~F' n B T . n ^G f n B t _ n _ l -H' n = o) 


SOLUTION OF FINITE-DIFFERENCE EQUATIONS 

Equations (10) and (11) represent two equations from 
which cr T , n and <x t , n can be expressed in terms of a T ,n-\ and 
<r t ,n-i ■ If the linear nature of the equations and the possi- 
bility of successive application of the equations to proceed 
from one station to the next are considered, the stresses at 
any station can ultimately be expressed in linear terms of 
the stresses at any other station. It will be convenient to 
express the stresses at all stations in terms of the stresses at 
the station a. At this station, the unknown value is the 
tangential stress <r t . a \ hence, the stresses at station n are 
expressed in the linear terms 

G T,n = A Ttn <T B t,ti 

G t,n^ = A t, n (T t, a -\- B t,n (12a) 

and those at station n— 1 in the form 

O’r.n— 1 == A r ,n—lG t,aA~ B r ,n — 1 

°'c»-l = -At p n- 10 'f, a -|- 2 ?«,n-l (12b) 

where the coefficients A r , n , B r , n , A t , n , and B ln are as yet to 
be determined. 

The substitution of equations (12a) and (12b) into equa- 
tions (10) and (11) and the separation of the terms with and 
without dt.a result in the equations 

(G n A Tn T) n At, n F n A rn _ i G„A t.n-l ) G t,a~\~ 


from which A r , n , A t ,„ } B r ,„, and B lH can be determined in 
the form 


A r n — K.nA r , n _\ -)- L n At, 

A t , n — K n A r , n —i~ L n A t, ti—i 
B r , n—K n B r , n _ x -\- L n B t , 

Bt, n — K' n B T , E' n D t, n-iFM', 


(16) 


If the coefficients A r ,„, A t , n , B r ,, n and B t , n are known for 
station n- 1, they can be determined by means of equation 
(16) for station n. 

The coefficients at the first station (r — a) can be deter- 
mined by inspection for both the solid disk and the disk 
with a central hole. Inspection of equation (12a) shows that 
for a solid disk in which both the tangential and radial 
stresses at the first station are equal to a ( , a 

A T , a — A t, a — 1 
Br, a =B t , a = 0 


For a disk with a central hole in which the radial stress 
at the first, station is zero and the tangential stress is <r t , a 


A t , a — Br, a B t, a 6 


A,. a = 1 


( C n B T , n T) n B t,n F n B T , n _ i — G n B H n ) — 0 (13) 

and 

( C'nA T , n — D'nA t ,n — F' n A T , n _iFG'n A /.„_]) G t , a -\- 
(C' n B T , n -D' n B t ,n-F' n B M +G' n B t , n _ l -H' n ) = 0 

(14) 

The stress a t , a is really arbitrary as far as equations (13) 
and (14) are concerned because it depends upon the boundary 
conditions and not on the equations of elasticity from which 
equations (13) and (14) were derived; that is, by a suitable 
choice of the factors that determine boundary conditions, 
such as blade loading and shrink fit., a t , a can be set at any 
desired value without invalidating in any way the equations 


From these known coefficients at the first station, the coeffi- 
cients at all other stations can be determined by successive 
applications of equation (16). Once all the coefficients have 
been determined, the unknown <r ta can be determined. 
The radial stress at the rim <r r , b is the centrifugal loading 
of the blades 

G’r, b — - A r , <J (, a “f~ Br,„ 


or 


G t. a * 


G r, b $t , 

A t , b 


h 


(17) 


where A r , b and B r ,i, are the coefficients for radial stress at 
the rim. The radial and tangential stresses at all stations 
can be obtained from equation (12a) after a ( , a and all the 
coefficients have been determined. 
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ILLUSTRATIVE APPLICATIONS 

CASE 1— ELASTIC-STRESS DISTRIBUTION IN SOLID DISK 

' The profile of a disk that is to be analyzed for stress dis- 
tribution at a speed of 11,500 rpm and the temperature 


I 



(a) Disk profile: Case I, solid disk; Case II, disk with 1-inch central hole. 

(6) Variation of temperature along disk radius. 

(c) Variation of Poisson’s ratio along disk radius. 

(rf) Variation of modulus of elasticity along disk radius. 

(e) Variation of coefficient of thermal expansion along disk radius. 

Kicukk 2.— Disk profile, temperature distribution, and variation of physical properties of 
disk material as function of radius, for illustrative problems. 


distribution of the disk are shown in figures 2 (a) and 2 (b), 
respectively. The first step in the analysis is to choose 
an arbitrary number of stations along the disk radius. The 
first station is chosen at a radius of about 5 percent of the 
rim radius, the last at the rim. The stations need not be 
equidistant; in fact, it is advisable to choose the stations 
closely together where there is sharp change in disk contour, 
in temperature gradient, or in variation of physical proper- 
ties. In this case 18 stations were chosen, spaced relatively 


close together near the rim where the gradients in tempera- 
ture and physical properties were high and near the center 
for subsequent use of the same example to illustrate the 
effect of a central hole. When only a solid disk is considered, 
no concentration of points near the center is necessary. 
The various steps of the calculation are tabulated in table I. 

The disk radius at each station is listed in column 1 of 
table I. The thickness of the disk at each station is listed in 
column 2. A sharp discontinuity in thickness, such as an 
abrupt flange, should be faired in the disk contour and the 
faired disk used in determining thickness. 

Ordinarily the density of the material is constant through- 
out the disk, even over the wide range of temperatures. If 
a faired disk has been used, however, the density of the 
material in the faired region should be adjusted to produce 
the total mass that actually exists in the region of each 
station. Although a flange does not reduce the stress at its 
own region by increasing the area, its mass must be included 
as it produces centrifugal stresses throughout the disk. The 
corrected density at each station multiplied by the squat e 
of the rotational speed is listed in column 3. In this case 
no fairing was necessary; hence, all values of density are 
equal. 

Poisson’s ratio, listed in column 4, has only an insignificant 
effect on the stress distribution and, because no accurate 
data are available, a constant value of 0.333 may be used. 
If accurate data on the variation with temperature of Pois- 
son’s ratio are available, use of the exact variable values 
presents no greater difficulty than use of a constant value. 
The values of ju used in this example are shown in figure 2 (c) , 
and were for convenience obtained by the assumption of a 
linear variation in m with temperature. 

The modulus of elasticity at each station is listed in column 
5. Variations in this property have a significant effect on 
the final stress values and accurate data should be used if 
available. For this example, E was arbitrarily assumed to 
depend linearly upon the temperature, and the variation 
along the radius is shown in figure 2 (d). In practical 
computations, the true values of elastic modulus associated 
with the particular temperature at each station may be used. 

The coefficients of thermal expansion are tabulated in 
column 6. These coefficients must be the average values 
applicable to the range between the temperatures actually 
existing and those at which there is no thermal stress. lor 
a homogeneous disk in which there is no shrink fitting of one 
part to another, the condition of zero thermal stress is at 
room temperature. Engineering tables usually list the 
average temperature coefficient of expansion between room 
temperature and values of high temperature, the listed \ allies 
may therefore be used directly. 

The difference between the actual temperature and the 
temperature at which there is no thermal stress is listed in 
column 7. In this case, the stress-free condition is at a room 
temperature of 70° F. Column 7 is therefore obtained by 
subtracting 70° F from each of the values in figure 2 (b). 
This column is of great significance in the case involving 
shrink fits. 
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TABLE I.— CALCULATION OF STRESSES IN DISK WITH NO CENTRAL HOLE 

[Engine speed, 11,500 rpm; operating temperature, 1270° F at rim, 670° F at center; average radial stress at rim, a r .b, 8500 lb/sq in.] 
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The quantities C n to M' n are computed for eaeh station as 
indicated in columns 8 to 34 of table I. Values in each of 
these columns can be obtained in one set of operations on a 
standard computing machine. The method of obtaining 
the data from the suitable previous columns is indicated at 
the heading of each column. 

Values in columns 33 and 34 must be simultaneously 
computed. The first value for each of these columns is 
unity. Subsequent values ' make use of the previously 
obtained values in the same columns. Thus, to determine 
the value for column 33 at station 2, column 27 at station 2 
is multiplied by column 33 at station 1, and the product is 
added to the product of column 28 at station 2 by column 34 
at station 1. For example: 

0.8 1 902 X 1 .0 +0. 18098 X 1 .0 = 1 .00000 

Columns 35 and 36 are likewise simultaneously computed. 
The first value in each of these columns is zero and each 
subsequent value is obtained from the previous values in 
accordance with the symbolic notation given at the head of 
each column. Thus, to obtain the column 35 at station 2, 
column 27 at station 2 is multiplied by column 35 iit station 1, 
column 28 at station 2 is multiplied by column 36 at station 1, 
and the two products arc then added to column 31 at 
station 2 


0.81902X0 + 0.18098X0 — 67.192 = — 67.192 

Column 37 is uniform for all stations and is obtained from 
the expression 


< ?t ,h — (35) ft 
(33) b 

where <r rJ> is the blade loading at the rim. The blade loading 
is obtained by dividing the total centrifugal force at the root 
of the blades by the total rim peripheral area. In this 
problem a r _„ is 8500 pounds per square inch; column 37 is 
therefore, 

8500 ,896 ) _ 42,008 pounds per square inch 

Columns 38 and 39 give the radial and tangential stresses, 
respectively, at each of the stations. As indicated in table I, 
they are obtained by routine multiplications and additions 
of columns 33 to 37. 

The radial and tangential stresses from columns 38 and 39 
are plotted in figure 3. The stresses at the center of the disk 
are taken equal to those at station a, which is % inch removed 
from the center. 

Because' the method presented is the only one known to 
the author that tak(fs into account point-to-point variation 
in Poisson’s ratio, the error involved in the assumption of a 
.constant value of this quantity as compared with the rigorous 
treatment of its point-to-point variations is valuable to 
determine. The broken-line curves in figure 3 show calcula- 
tions for constant values of /*=0.3 and /x = 0.5 compared 
with the solid curves, which show the stresses for a contin- 



Fkjure 3.— St.ross distribution in solid disk of figure 2 computed with constant and with 
variable values of Poisson’s ratio. 


uously variable value of n with temperature, as shown in 
figure 2 and tabulated in column 4 of table I. The near 
coincidence of these curves indicates that the assumption of 
a constant value of n within the range of actual values results 
in accurate final values of radial and tangential stresses. 

The effect that difference in the number of stations has 
on the accuracy of the results is shown in figure 4. Little 
accuracy is gained by the use of additional points; as few 
as six points in this particular case can yield accurate results 
at a great saving in computing time. 

CASE II— ELASTIC-STRESS DISTRIBUTION IN DISK WITH CENTRAL HOLE 

A disk with a central hole is studied in a manner similar to 
the solid disk except that the first station is taken at the in- 
side boundary instead of at an arbitrary small distance as in 
the solid disk. The choice of stations near the central hole 
is, however, critical for this case. Stations should be taken 
also at distances of 1, 2, 3, and 5 percent of the rim diameter 
from the inside boundary of the disk. In order to illustrate 
the procedure, the disk of figure 1 is again used but with a 
central hole 1 inch in diameter, chosen so that station a will 
be conveniently located in the same place as station a 
for Case I. 

Columns 1 to 32 for the disk of figure 2 (a) with the central 
hole are identical to the corresponding columns of table I 
for the solid disk. In column 33 the entry for station a is 0 
instead of 1 as for the solid disk; otherwise, the procedure for 
calculating columns 33 to 39 is the same as that of table I. 
Table II gives the modified columns 33 to 39 that result from 
changing the single first entry in column 33 and figure 5 
shows a plot of the resulting radial and tangential stresses; 
the curves marked “18 stations” are the stress values for 
this computation. 
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Figure 4.— Effect on calculated stress distribution in disk of figure 2 of various numbers of 

stations. 

The results of supplementary calculations using different 
numbers of stations (fig. 5) indicate that considerable error 
can result in the determination of the peak stress at the inner 
boundary if an insufficient number of. stations are chosen 
near this boundary. A more judicious choice of stations for 
the 6- and 10-station systems could produce more accurate 
results than those shown; in the absence of experience in 
choice of locations, however, it is better to choose a large 
number of stations and insure accuracy. 

A practical procedure used to reduce the amount of cal- 
culation necessary to obtain the critical end stresses is to 
calculate the stress distribution on the basis of a solid disk 
using a few stations and then to modify the stresses in 
the immediate vicinity of the central hole by the stress- 
concentration factor characteristically introduced by the hole. 
A comparison of figures 4 and 5 indicates that, with the ex- 
ception of the region immediately adjacent to the central 
hole, the stresses are similar for the cases of the solid disk 
and of the disk with the central hole; this stress distribution 
depends very little on the number of stations chosen. By 
reference II (fig. 145), for example, the characteristic stress 
concentration for a disk with a central hole of which the 
diameter is one-twentieth of the outside diameter is about 
2.0. From calculations based on different numbers of 
stations, the calculated average stress at the center of the 


solid disk (fig. 4) is 43,000 pounds per square inch. The 
tangential stress at the inside boundary for the disk with the 
central hole should therefore be 2X43,000=86,000 pounds 
per square inch. The radial stress at a free boundary is, 

TABLE II— CALCULATION OF STRESSES IN DISK WITH 
1-INCH CENTRAL HOLE 


[Engine speed, 11,500 rpm; operating temperature, 1270° E at rim. 670° F at center] 
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Figure 5.— -Effect on calculated stress distribution in dis.k of figure 2 with central hole of 
various numbers of stations. (Location of stations shown on fig. 1.) 
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of course, zero. A curve faired between these boundary 
values and the general curves of figure 4 would coincide very 
elosely with the 27-station result of figure 5. 

CASE III— ELASTIC-STRESS DISTRIBUTION IN COMPOSITE WELDED DISKS 

For some applications, turbine wheels must be fabricated 
by welding parts composed of several materials. The 
method of analysis presented is applicable to studies of 
composite welded disks in which various alternatives of 
boundary location and shrink interference can be investi- 
gated with few changes in the tabulated computations. The 
procedure is illustrated for a typical application in which 
the boundary location is constant. 

The disk of figure 2 (a) is assumed to be made in two parts 
with the boundary at the 6-inch station. Figure 6 shows 
the two portions of the disk just before welding. The 
heat-resisting outer portion is heated to 670° F while the 
inner portion is maintained at 70° F. (In practice both 
portions may be heated while maintaining a desired tempera- 
ture differential.) At this temperature condition, an exact 
fit exists between the mating tips of the two parts. The 
wedge is then filled with weld metal. 

The assumption is made that this temperature differential 
between the two portions of the disk is maintained throughout 
the welding process in making the calculations. Any cooling 
of the outer region prior to the placement of the weld metal 
would produce a crushing of the mating tips and reduce the 
effective amount of shrink. Localized effects of the weld 
metal in producing residual stresses are neglected. The 
calculations are made as if the high-temperature alloy, hav- 
ing full width at the mating face, is shrunk at 670° F onto 
the full-width steel central portion. Table III shows the 
essential tabulations for this case. 

In order to insure accuracy, a few more stations than were 
used in tables I and II have been chosen in the vicinity of 
the boundary. The densities of the two materials are some- 
what different. The quantities for ij. and E are the values 
at the temperatures of figure 2 (b). The quantity AT at 
each station is the difference between the existing tempera- 
ture and the temperature at which there is zero thermal 
stress. The temperatures of zero thermal stress occur just 
before the shrink fit when the outer portion is at 670° F and 
the inner portion is at 70° F. Therefore, for the outer 
portion 670° F is subtracted at each station from the tem- 
peratures of figure 2 (b) and for the inner portion 70° F is 
subtracted. The value of a at each station must be the 
average a between the stress-free temperature and the operat- 
ing temperature. At the rim., for example, an average 
coefficient of expansion between 670° and 1270° F must be 
used. The average temperature coefficient o> _ 2 applicable to 
the range between any two temperatures T x and T 2 can be 
found by 

ai - 2 ~(T 2 -TMl + a\T\) 

where' a', and a'-, are*, the* average' coefficients of thermal ex- 
pansion between room temperature and the temperature's 
7’, and T>, respectively, and 7 T/ , and 7L are the temperature 
differences between 1\ and T> and room temperature. 


The procedure of calculating table III from column 8 on 
is similar to that of table I. The final calculated values of 
stress are shown in figure 6. 

Comparison of figures 5 and 6 shows that a shrink fit, 
unless excessive, can have beneficial effects. The shrink fit 
reduces the tangential tensile stresses that exist near the 
central hole under operating conditions and also the tangen- 
tial compressive stress at the rim. Compressive stresses at 
the rim can be detrimental. If the elastic compressive stresses 
exceed the yield strength of the material, plastic flow takes 
place and a residual tangential tensile stress exists after 




tvai suction, 70° K. 

operation. Because the region of the rim is a stress- 
concentrated area as a result of the blade attachments, even 
relatively small residual tensile stress may cause cracks. 
The shrink fit removes the high tensile stress at the center 
and the high compressive stress at the rim but introduces a 
high tensile stress at the boundary of the two fitted regions. 
The boundary is at a lower running temperature and has no 
stress-concentrating effects of the blade attachments. 

The optimum amount of shrink, however, is fairly critical. 
Probably the shrink of the illustrative example is excessive. 



TABLE III— CALCULATION OF STRESSES IN SHRUNK DISK 

[Engine speed, 1 1,500 rpm; operating temperature, 1270° F at rim, 670° F at center; shrinking condition, 670° F for rim portion, 70° F for central portion] 
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There is no need to reduce the stresses at the hole and at the 
rim as much as shown in figure 6 at the expense of such a 
high stress at the boundary. An additional calculation can 
readily be made by using a smaller temperature difference 
of shrinking than assumed in this calculation. Only col- 
umns 6, 7, 24, 31, 32, and 35 to 39 are affected by any change, 
and the redistribution of stress can be calculated very, rap- 
idly. Thus, a more suitable shrink fit can readily be found. 

CASE IV— CHECK ON ADEQUACY OF METHOD 

Checks on the adequacy of the method were obtained by 
comparing the results of finite-difference calculations to 
theoretically correct results in several cases where the latter 
could be obtained. In one case a parallel-sided disk was 
studied. The conditions of operation are shown in figure 7; 


120 x 10* 



Figure 7.— Comparison between theoretical and finite-difference-solution stresses in parallel- 
sided disk of 20-inch diameter rotating at 10,000 rpm with temperature gradient that varies 
as fourth power of radius from 600° F at center to 1200° F at rim. K, 30X10 4 pounds per 
square inch; a, 8X10~ 9 (in./(in.H°F)); n, 0.3. 


the circles and squares show the radial and tangential 
stresses as determined by the finite-difference method, 
respectively; and the solid lines show the theoretically 
correct stresses obtained by rigorous solution of equations 
(1) and (7) for this case in which E, a, and /x are constant. 
This correlation is seen to be very good. The maximum 
deviation occurs at the boundary of the central hole where the 
difference between the tangential stress as computed by the 
finite-difference method and the theoretically correct value 
is about 2 percent. The average deviation between the 
theoretical and finite-difference stresses throughout the disk 
is less than ±K percent. Checks on a solid disk produced 
closer agreement, even when a small number of stations was 
used. A check on a disk of uniform strength produced 
results differing from the exact solution in the order of ± / 
percent of the theoretical stresses throughout the entire disk. 

CONCLUSIONS 

The finite-difference method of calculating stresses in 
rotating disks has been applied extensively to various types 
of turbine disk under different conditions of constant tem- 
perature or with a temperature gradient. The procedure 
was found to be convenient and rapid. Where checks were 
available, the results showed a high degree of accuracy. 


Flight Propulsion Research Laboratory, 
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